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Universal Fermi Gas with Two- and Three-Body Resonances 
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We consider a Fermi gas witli two components of different masses, with the s-wave two-body 
interaction tuned to unitarity. In the range of mass ratio 8.62 < M/m < 13.6, it is possible for 
a short-range interaction between heavy fermions to produce a resonance in a three-body channel. 
The resulting system is scale invariant and has universal properties, and is very strongly interacting. 
When M/m is slightly above the lower limit 8.62, the ground state energy of a 2:1 mixture of heavy 
and light fermions is less than 2% of the energy of a noninteracting gas with the same number 
densities. We derive exact relationships between the pressures of the unitary Fermi gases with and 
without three-body resonance when the mass ratio is close to the critical values of 8.62 and 13.6. 
Possible experimental realization with cold atoms in optical lattices is discussed. 

PACS numbers: 05.30.Fk, 03.75.Ss 



Introduction. — Dilute Fermi gas at infinite scattering 
length P, 0, Q has attracted considerable attention re- 
cently. The system can be realized in atomic traps using 
the Feshbach resonance 0, H, 0, H, This system 
is called Fermi gas at unitarity and has interesting uni- 
versal properties. For example, the energy of the ground 
state of a symmetric Fermi gas at unitarity is £^ = f -Efroc, 
where i?free is the energy of a noninteracting system with 
the same density and ^ is a dimensionless number. Most 
recent evaluations for ^ yield ^ « 0.4. 

Other systems with universal properties that have been 
considered so far are variations of this unitary Fermi 
gas. One can consider different densities of the two 
fermion components lO^llL or different masses of the 
two fermions, or both 1^, ll3|. Here we show that there 



exists a novel type of system with both two- and three- 
body interactions fine-tuned to resonance. 

Three-body wave function. — In order to realize a three- 
body resonance, it is necessary to allow the two compo- 
nents of the Fermi gas to have different masses. We shall 
denote the mass of the heavy fermion as M and of the 
fight fermion as m. The ratio u = M/m is a free param- 
eter. 

Let us recall what is known for three-body physics, 
using the Born-Oppenheimer approximation at first [3] , 
but then proceed to the exact result. Consider a three- 
body system of two heavy and one fight fermions. If M ^ 
771 one can use the adiabatic approximation, separating 
the wave function into the fast part of the light fermion 
and the slow part of the heavy fermions, 

%l:{xi,X2,y) = Vsiow(iCi,a;2)7/'fast(y;a;i,a;2). (1) 

The function ';/'tast(y; a^i, 3^2) satisfies the Schrodinger 
equation of one particle moving in the field of two static 
centers: 

- 7;—^'ii}fast{y\Xi,X2) = E{Xi,X2)^/j{a.st{y)- (2) 

The resonant two-body interaction implies that the 
asymptotics of 7Afast(y) when y Xi, i = 1,2 is 



\y — Xi\ ^ -\- 0{\y -~ Xi\). One finds one negative energy 
level with 
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, R=\xi-X2\ 



(3) 



where C — 0.5671 ... is the solution to equation e^^ — x. 
We then regard E{R) as the potential energy between the 
two heavy fermions. Since it is a potential, the wave 
function ■0siow near the origin behaves as R^YimiO, cj)), 
where a is a solution to the equation 



a{a I) = 1(1 + I) - 



(4) 



Because of the Fermi statistics I is odd. When the right- 
hand side of Eq. (jH) is greater than —j, there are two real 
solutions a+ and q;_ , relating to each other by + «_ — 

— I (a+ > — 1/2 > «_). If the right-hand side is less than 

— J, which is the case for Z = 1 and u > Mmax ~ 14.0, 
then a is complex. This corresponds to the Efimov effect, 
where the system develops deep bound states in the zero- 
range limit, which break scale invariance. 

In the former case, u < Wmax, there are two further 
possibilities. If a+ > i and a„ < — f , then the solution 
V'siow R°'~ is not normalizable at the origin. There- 
fore, the only consistent boundary condition for V'siow 
at the origin is tAsIow However, when < i, 

Q!_ > — |, both asymptotics are normalizable at the 
origin. Therefore, generically ipsXow will behave as 

7^slow(i?,^,</') - (c+i?"+ +C_i?"-)rz„(0,(/)), (5) 

where the ratio Cj^/c-^ is determined by the interaction. 
For a generic interaction c+ /c_ , which has dimension 



[length]"- will be of order 



where tq is the 



characteristic interaction range. That means that at dis- 
tances R^ ro the wave function behaves like . How- 
ever, by including an attractive force between the heavy 
fermions (in addition to the effective attraction induced 
by the light fermion), and fine-tuning its strength, one 
can achieve the situation where only the coefficient of the 



2 



R°'- piece is very large. This corresponds to a three-body 
resonance, and c+ = corresponds to this resonance be- 
ing at threshold. For ^ = 1, this is possible when u > 
Umin ~ 7.77. To verify our point, we considered a model 
potential V{R) = (fi^/M) {^i/R^'^ - (32/R^) between 
the heavy fermions. Within the Born-Oppenheimer ap- 
proximation, when P2/Pi^^ ~ 6.5 + 2.9a+ (— i < a+ < 
^), which is not a strong enough attraction to produce a 
two-body resonance, the desired three-body resonance is 
achieved. 

The Born-Oppenheimer approximation employed 
above is not entirely justified, because the obtained val- 
ues of Umin and Mmax are not parametrically large. The 
conclusions, however, are confirmed by an exact treat- 
ment of the three-body problem. The wave function 
'if{xi,X2,y) satisfies the free Schrodinger equation with 
the two-body boundary condition when |a;,; — y| — > 0. In 
hyperspherical coordinates, its behavior is 



(6) 



where R is the hyperradius and Q denotes the hyperan- 
gular variables, and 7 is related to a by 7 = a — |. For 
the I = 1 channel, the power 7 is determined from the 
equation |15|] 

cos[(7-|-l) arccos :;^] — cos[(7+l) arccos ■^^] 



(7 -I- 1) sin7r7 



+ 



'[(7+3) i 



ti+iJ 



'[(7+3); 



(7 -|- 3) sin7r7 



2u' 



rV2^m. (7) 



[u + iy 

Solving this equation, we find that w,nax ~ 13.6 and 
Umin ~ 8.62, which somewhat differ from the values ob- 
tained in the adiabatic approximation. At the qualita- 
tive level, however, the physics is the same: for the mass 
ratio between these two values, one can fine-tune one 
three-body channel to resonance. Note that the mass ra- 
tio Winin ~ 8.62 played a special role in the discussion of 
three-body recombination in Ref. llQ. 

Many-body physics. — We now discuss the ground state 
energy of the many-body system with both two- and 
three-body resonances. The first question is whether the 
system is stable toward collapse. Such a collapse can hap- 
pen if there exists an Efimov state of four or more parti- 
cles for the mass ratio between 8.62 and 13.6. This would 
be very unusual, as it is before the three-body Efimov 
state appears. Consider four-body states. In the Born- 
Oppenheimer approximation, there is no Efimov state of 
two heavy and two light fermions, since there is only one 
bound state in the Schrodinger equation with two static 
centers. The case of three heavy and one light fermion 
is more subtle, but here the three centrifugal potentials 
between the heavy fermions make the Efimov effect less 
likely. We shall proceed under the assumption that there 



is no Efimov effect for any number of heavy and light 
particles, and the system is stable. 

Let us define ^ as the ratio of the ground state energy 
E to the ground state energy Ef^ee of the noninteracting 
system at the same densities of the heavy and the light 
components, rih and n;, respectively. 



E{nh, ni) = ^Efrccinh^ni) 



(8) 



Since the system has no intrinsic dimensionful parameter, 
^ can depend only on two variables: the mass ratio M/m, 
and the relative abundance of the two components. For 
the latter, we define 



ni 



rih + ni 



(9) 



In general, the problem of finding ^(M/m, a;) is diffi- 
cult, since there is no small parameter in the problem. 
We shall concentrate on the case when M/m is slightly 
larger than the lower critical value 8.62. In this case, the 
normalization integral for the three-body wave function 
is dominated at small distances. Thus, one can think 
about the three-body system as a localized bound state. 
The situation is analogous to the behavior of the two- 
body state at resonance near four spatial dimensions. As 
the number of spatial dimensions d approaches four from 
below, the two-body wave function becomes more and 
more dominated by short distances, which has led Nussi- 
nov and Nussinov to suggest that ^ ^ in this limi t Il6l . 
This was confirmed by the e expansion near d = i [l7[ . 

The view that the three-body system behaves as a 
tightly bound state is confirmed by its behavior in an 
isotropic harmonic potential. It is known that each 
exponent 7, determined from solving the free-space 
Schrodinger equation, is related to an energy level of the 
three particles in an isotropic harmonic potential [l8| . 



Eo^ ( ^ + 7 ) ^w, 



(10) 



where cu is the oscillator frequency. When 7^—3, the 
energy approaches ^huj. This is consistent with the pic- 
ture of the trimer being a localized bound state of all 
three particles. We emphasize that the three-body state 
is not a real bound state with a finite binding energy. 
In particular, in the harmonic trap there are breathing 
modes, corresponding to a ladder of energy excitations 
Eq + 2nhw, n = 1, 2, . . . [loj, that are absent for the sin- 
gle particle. However, when the trimer is in the ground 
state with respect to the breathing mode, it behaves as 
a pointlike particle. 

Having in mind this picture of trimers as objects local- 
ized in space, let us put an upper bound on the value of 
^ at X = ^, and u = 8.62 -I- e. In this case one can com- 
bine all heavy and light particles (whose relative number 
is 2:1) into trimers, and construct a dilute gas of such 
trimers. This gas will be a weakly interacting Fermi gas 
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of trimers with three possible polarizations (the trimer 
has 1 = 1). Its energy density is therefore 



5/3 



= c (11) 

V 32/3(2M + m)' ^ ^ 

where c = ^(GTr^)^/-^?!^ and ni is the density of the light 
fermions. On the other hand, if one turns off the inter- 
action, the system becomes a Fermi gas of light fermions 
with density rii plus a Fermi gas of heavy fermions with 
density = 2ni. The energy density of this gas is 



Efiee 

c — ' — 

V m 

Taking the ratio we find 



M 



1.93 X 10" 



(12) 



(13) 



Efrcc 32/3(u + 25/3)(2m+1) 

It is clear that the trimer gas energy is not the lowest 
ground state energy and hence is only an upper bound 
on the latter. For example, a trimer with energy near the 
Fermi energy can decay into its constituents (two heavy 
and one light fermions) since the binding energy is zero, 
and this decay can happen until chemical equilibrium 
is reached. The true ground state is therefore a mix- 
ture of the trimer gas with a gas of constituent fermions. 
Because of the strong interactions between constituent 
fermions, we do not have an exact formula for the ground 
state energy. The upper bound, on the other hand, im- 
plies that the system is very strongly interacting: the 
ground state energy is less than 2% of the energy of the 
noninteracting system: ^{u ^ Umin, i) < 1.93 X 10~2. 
This can be contrasted with the case of symmetric uni- 
tary Fermi gas with equal masses, where ^ « 0.4. 

For a general value of x, we can put an upper bound on 
^ by considering the state where the maximal number of 
trimers is bound and form a trimer Fermi gas, while the 
remainder (heavy fermions for x < ^ and light fermions 
for a; > i) form another free Fermi gas. The resulting 
bound is plotted in Fig. [1] One sees that the variational 
bound is lowest when x is close to 1/3. From the point of 
view of the ground state energy, the 2:1 mixture of heavy 
and light fermions is therefore the most interesting. 

Exact relationships near u = 8.62 and 13. 6. — For mass 
ratios u between 8.62 and 13.6, therefore, there exist two 
scale invariant regimes, which give rise to unitary Fermi 
gases with and without three-body resonance. We now 
show that near the lower critical value u^un there exists 
an exact relationship between the pressures of a unitary 
Fermi gas with three-body resonance, P, and a unitary 
Fermi gas without three-body resonance, Pq, at the same 
chemical potentials of light and heavy fermions, fih and 
fj,i. Namely, 



lim [P{ni,fih) - Po{ni,fih)] 

— >Umin 

_ [2(2M-t-TO)]3/2 



(2Ai„ + M/)'/'. (14) 



1.0 r 

0.8 - 

0.6 - 

0.4 - 

0.2 - 



0.0 0.2 0.4 0.6 



0.8 



1.0 



FIG. 1: The variational bound on ^ as a function of the abun- 
dance of light fermions x at M/m = 8.62 + e. 



We emphasize here that we have in mind the limit when 
for each mass ratio u we take the zero-range limit, com- 
pute the pressures P and Pq, and then take the limit 
u —^ Wmin- In general the limit u — s- Umin and the zero- 
range limit do not commute. 

The right-hand side of Eq. is the pressure of an 
ideal Fermi gas of trimers, with Z = 1, at chemical poten- 
tial 2fifi+fM- Physically, Eq. (jl4p means that at m — > Umin 
the unitary Fermi gas with three-body resonance con- 
sists of two parts: a trimer gas and a Fermi gas of the 
fermionic constituents outside the trimers. These two 
parts interact weakly with each other, but are in chemi- 
cal equilibrium. 

When a; — i, there is another simple relationship be- 
tween the ground state energies with and without three- 
body resonance. ^ of the Fermi gas with three-body reso- 
nance is determined by the analogously defined quantity 
of the unitary Fermi gas without three-body resonance 
Co via 



-3/2 



Co 



-3/2 



374. {x 



d- (15) 



We have argued above that the trimers decouple at 
Wmin- What is less trivial is that the fermions outside 
the trimers contribute the same pressure as a Fermi gas 
without three-body resonance. This seems strange: the 
three-body wave functions satisfy different boundary con- 
ditions at small distances {R^^ vs. R''+) in the two cases. 
To see that it is the case, consider the three-body wave 
function of three particles at some positive energy. The 
radial part of the Schrodinger equation reads 



^^^p 5 dij -f{-f + 4) 

+ R'dR ^2 

This equation has two solutions: 



^ = _fc2^(i?). (16) 



(17) 



where J-i, corresponds to the case with three-body res- 
onance, and Jjy to that without the resonance. The limit 
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tl/th 


8.62 


9 


10 


11 


12 


13 


13.6 


W/tn 


-3.80 
±0.03 


-3.625 
±0.025 


-3.15 
±0.05 


-2.575 
±0.075 


-2.0 
±0.1 


-1.2 
±0.1 


> -0.4 



TABLE I: Approximate values of W/th for three-body reso- 
nance on a lattice. 



u — > Mmin corresponds to v ^ 1. By using the relation- 
ship between the Bessel functions, 

Jv{z) cos{vk) — J^u{z) = sin(i^7r)Y'^(z), (18) 

one sees that in the limit v ^ 1 the two wave func- 
tions in Eq. (|17p are the same up to a sign, except for 
a very small region near the origin, kR - ^/^^. This 
means that, for the fermions outside the trimers, which 
have finite kinetic energy, the resonant and nonresonant 
boundary conditions are equivalent. This leads to the 
relationships (|14p and (fT5| . 

On the other hand, the limit u u^aa^ corresponds to 
i/ ^ 0, where two solutions in Eq. (|17p coincide. There- 
fore there is no distinction between Fermi gases with and 
without three-body resonance and we obtain another ex- 
act relationship 



lim [P{^j,i,fj.h) - PoinhUh)] = 0. 



(19) 



Other channels. — So far we have considered only the 
I = 1 three-body channel. One may ask whether similar 
fine-tuning can be done in other channels. To answer this 
question, one has to compute the exponents 7+ in other 
channels and look for the mass ratio u where it crosses 
— 1. This crossing happens only for odd I. For / = 3 
it occurs at u — 70.1, and the value of the mass ratio 
increases with increasing I. For these large mass ratios, 
the Efimov effect has already taken hold in the I — 1 
channel, and the system cannot be scale- invariant. We 
conclude that scale invariance with resonant three-body 
interaction can be achieved only in the I = 1 channel. 

Experimental realization in a lattice. — Fermions of 
species c and 6 on a cubic lattice with the Hamiltonian 



<ij> i 

<ij> 

(where extends over nearest- neighbor lattice sites) 

approach the above universal Fermi gas limit when the 
average number of fermions per lattice site <C 1, if param- 
eters are tuned to get cb and cch resonances: U/{th+ti) « 
—3.957, 8.62 <ti/th < 13.6, and numerically determined 
W is shown in Tab. U (currently we only know a lower 
bound for W at ti/th = 13.6). One can easily show that 
p-wave cc resonance is not reached until W/th ~ —9.53. 

The model (j20p may be experimentally realized with 
cold atoms on an optical lattice. If the cb scattering 



length in free space is small and negative, we can tune 
cb to resonance by tuning the lattice depth U/th 
is also tunable. Nearest-neighbor attraction between c 
atoms may be mediated by a third kind of atom having a 
much smaller effective mass than c and b, interacting with 
c atoms with an effective scattering length Ocs: < 
and jocal slightly exceeds the lattice spacing. Tuning 
with a magnetic field near a c3 Feshbach resonance, we 
can tune the magnitude of W in Eq. (PO)) . The lattice 
also dramatically increases the lifetime of the system by 
suppressing the three-body recombination rate 21 1 . 

Conclusion. — We have argued that there exists a class 
of universal Fermi gases with resonances in both two- 
and three-body channels. This system remains strongly 
coupled in the whole range of mass ratio, 8.62 < M/m < 
13.6, where it exists. We have discussed the possibility of 
realizing these gases with cold atoms in an optical lattice. 
The relevant lattice model, Eq. (PO)) . may also be studied 
with quantum Monte-Carlo simulations. 
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